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We present new measurements of the 7 Be(p,7) 8 B cross section from E cul = 116 to 2460 keV, that 
incorporate several improvements over our previously published experiment, also discussed here. 
Our new measurements lead to Si7(0) = 22.1 ± 0.6(expt) ± 0.6(theor) eV b based on data from 
E crn = 116 to 362 keV, where the central value is based on the theory of Descouvemont and Baye. 
The theoretical error estimate is based on the fit of 12 different theories to our low energy data. 
We compare our results to other Si7(0) values extracted from both direct ( 7 Be(p,7) 8 B) and indirect 
(Coulomb dissociation and heavy-ion reaction) measurements, and show that the results of these 3 
types of experiments are not mutually compatible. We recommend a "best" value, Si7(0) = 21.4 ± 
0.5(expt) ± 0.6(theor) eV b, based on the mean of all modern direct measurements below the 1 + 
resonance. We also present S-factors at 20 keV which is near the center of the Gamow window: the 
result of our measurements is Si7(20) = 21.4 ± 0.6(expt) ± 0.6(theor) eV b, and the recommended 
value is S ir (20) = 20.7 ± 0.5(expt) ± 0.6(theor) eV b. 

PACS numbers: 26.20+f, 26.65+t, 25.40Lw 



I. INTRODUCTION 

Studies of high-energy neutrinos produced by 8 B de- 
cay in the sun have shown that when these neutrinos are 
detected on earth the number that remain electron neu- 
trinos (^ e 's) is only about half of the number that have 
oscillated into other active species |l| . Precise predictions 
of the production rate of 8 B solar neutrinos are important 
for testing solar models, and for limiting the allowed neu- 
trino mixing parameters including possible contributions 
of sterile species. 

The predicted 8 B production rate Q is based on solar- 
model calculations that incorporate measured reaction 
rates for each of the solar burning steps following the 
initial p + p reaction, the most uncertain of which is 
the Be(p,7) 8 B rate. The currently recommended value 
for Si7(0), the astrophysical S-factor for this reaction, is 
19 ±2 eV b Q. However, Si7(0) must be known to better 
than ±5% in order that its uncertainty not be important 
in the predicted v e production rate [2J . 

Because of the key role of Si7(0), it has been measured 
many times using a variety of techniques: direct stud- 
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ies of the 7 Be{v^A s B cross section S IE 0, H H E3 
ITU IT2I lip. Ibl flal. and indirect studies using Coulomb- 
breakup |l6lll7lll8l|l9l.l20f or peripheral heavy-ion trans- 
fer and breakup [2lL l22j| reactions. The direct technique 
has the advantage that it studies the reaction that ac- 
tually occurs in the sun, but is difficult experimentally 
because it requires a radioactive target or beam, and the 
cross sections are small. 

Our direct measurements reported here incorporated 
several improvements over traditional methods. We elim- 
inated a major uncertainty in many previous experiments 
due to uncertain and nonuniform target areal density by 
using a ~1 mm-diameter beam magnetically rastered to 
produce a nearly uniform flux over a small ~3.5 mm di- 
ameter target. We directly measured the energy-loss pro- 
file of the target using a narrow 7 Be(a,7) n C resonance. 
We made frequent in situ measurements of the 7 Be tar- 
get activity to determine the target sputtering losses. We 
also made the first (and, to date, the only) measurements 
of 8 B backscattering losses. 

This paper reports the results of three separate exper- 
iments that determine the 7 Be(p,7) 8 B cross section over 
a range of 3 orders of magnitude. The first 01 use d a 
target of 106 mCi initial activity (here called BE1 [2]|) to 
measure the 7 Be(p,7) 8 B cross section at mean center-of- 
mass proton energies, E cm = 186 to 1203 keV, including 
the Ml resonance near E cm = 630 keV. The a-detector 
solid angle in this experiment was determined using a 
7 Li(d,p) 8 Li reaction yield ratio. In a second, abbrevi- 
ated measurement with a new 7 Be target of similar activ- 
ity (called BE2), we determined the cross section in the 
range E cm = 876 to 2459 keV. In this second measure- 
ment, we became concerned over possible inaccuracies in 
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the 7 Li(d,p) 8 Li reaction yield ratio method for determin- 
ing the a-detector solid angle, which relied on calculated 
extrapolations of the continuous a-spectrum tail below 
the experimental cutoff (a similar but smaller extrapo- 
lation is necessary to interpret the 7 Be(p,7) 8 B reaction 
yields). 

Recently we completed a 7 Be(p,7) 8 B measurement 
with a 340 mCi target (called BE3) that incorporated 
several improvements over the BE1 experiment, includ- 
ing a thinner Si surface-barrier detector for which the cor- 
rection (and hence the uncertainty) due to the fraction 
of the 8 B /3-delayed a-spectrum lost below the detector 
threshold was minimal. We avoided using the 7 Li(d,p) 8 Li 
reaction for solid angle normalization, using instead a 
custom-made 1 Gd a-source fabricated on a backing of 
the same design as the one used for the 7 Be target. As 
a result, we were able to measure and minimize all im- 
portant sources of systematic error in determining the 
7 Be(p,7) 8 B cross section. This BE3 measurement, which 
covered the range E cm = 116 to 1754 keV, yields our 
best determination of the absolute cross section. Hence 
we base our absolute S-factor determination on this ex- 
periment. Our new determination of Si7(0) turns out to 
be in excellent agreement with our previously published 
value 0. 

We discuss all 3 experiments in Secs.lTTl- HVI below. with 
an emphasis on the BE3 measurements. Comments on 
the BE1 and BE2 experiments are labeled explicitly; un- 
labeled comments refer to the BE3 experiment. We esti- 
mate the extrapolation uncertainty in Si7(0) in Sec. IV Al 
and make detailed comparisons with other direct experi- 
ments in Sec. IV Bl We discuss indirect measurements in 
Sec. IV CI and in Sec. IVII we recommend a "best" value 
for Si7(0). We summarize our results in Sec. lVIll and we 
discuss their implications for solar and neutrino physics 
in Sec. IVTTT1 



II. EXPERIMENTAL PROCEDURE 

A. Target fabrication 

Our 7 Be targets were made at TRIUMF using a tech- 
nique described previously 0> H3 ■ Briefly, the 7 Be was 
produced by the 7 Li(p,n) 7 Be reaction at E p = 13 MeV, 
using a Li metal target and the TR13 cyclotron at TRI- 
UMF. The 7 Be was chemically separated using a glass 
filter and various solutions including HC1, then heated 
and dried, leaving the 7 Be predominantly in the form 
of 7 BeO. Then a 2-step reduction/distillation was per- 
formed in vacuum. In the first step, the 7 BeO was heated 
in a Zr-lined Mo crucible, depositing 7 Be metal on a Mo 
piece which served as the lid of the apparatus. In a second 
distillation the Mo piece served as the heating crucible 
and the 7 Be was deposited on a target backing. This 
technique yielded 'Be targets of high-purity (40 - 63% 
by atom number - see Sees. Ill G*l and IiTK|) . with initial 
activities of 106 mCi (BE1), 112 mCi (BE2) and 340 



mCi (BE3). The 7 Be activities produced at TRIUMF 
for these targets were roughly 220, 420 and 630 mCi, 
respectively. 

The target backings consisted of 1.3 by 1.5 cm Mo 
plates with stainless-steel water-cooling tubes brazed 
onto their back sides, and mounting brackets on the top 
for attaching the backing to the rotating arm (see Fig. 
4 of [25j). The front face of the plate consisted of a flat 
surface with a 4 mm diameter, 1.5 mm high Mo post in 
the center. Prior to 7 Be deposition, a Mo washer was 
press- fitted around the post. The washer and post were 
machined to precise tolerances to ensure a tight fit, and 
the post-plus-washer was machined flat after assembly. 
After 7 Be deposition, the washer was broken away, leav- 
ing 7 Be only on the post. This procedure ensured that 
7 Be remained only within a small «3.5 mm diameter cen- 
tral feature that was important for our large-area- 
beam/small- area-target technique described below. 

A 7-activity scan of the BE3 target is shown in Fig. ^ 
The scan was measured with a Ge detector collimated 
by a 51 mm thick "heavy metal" (tun gste n alloy) block 
containing a 0.125 mm slit (see also ref. |23). Scans mea- 
sured both before and after the 7 Be(p,7) 8 B experiment 
showed very similar shapes. Scans measured with the tar- 
get turned sideways showed that an insignificant amount 
of 7 Be was located elsewhere than on the top of the post, 
and were used to determine t he p osition-resolution func- 
tion for a line source (see ref. |25j). Fits to the measured 
scans using this resolution function showed that the Be 
density distribution on the post was constant within a 
radius of 1.5 ± 0.1 mm and decreased to zero at a radius 
of 1.8 ± 0.1 mm. 

Additional information on target purity and composi- 
tion is given in Sees. Ill Gl and III Kl below. 



B. Experimental apparatus 

The 7 Be(p,7) 8 B cross sections were measured at the 
University of Washington Center for Experimental Nu- 
clear Physics and Astrophysics, using p, d and a beams 
at energies up to 3 MeV from the tandem Van de Graaff 
accelerator with a terminal ion source. For the BE1 and 
BE2 experiments, proton beam currents up to 16 /iA were 
used. For the BE3 experiment a straight-field accelerat- 
ing tube was installed in place of the first high-energy 
spiral-field tube section to obtain lower energies and pro- 
ton beam currents up to 35 fiA. Typical a beam currents 
were 8 /it A or less on target. Based on thermal tests of 
our water-cooled backings, beam-power levels on target 
were kept to less than 10 watts to prevent target damage 
from beam heating. 

Figure |21 shows a top view of the target chamber. The 
beam, indicated by the shaded taper, enters the cham- 
ber from the upper left, and passes through a large-area 
aperture followed by a 31 cm-long cylindrical LN2 cold 
trap with inner and outer diameters of 2.9 and 11.4 cm. 
The cold trap had a removable Cu liner which captured 
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FIG. 1: BE3 Activity scan measured before (top panel) and 
after (bottom panel) the (p, 7) experiment. Solid curves - 
calculated activity (see text). 



most of the 7 Be sputtered from the target during long 
bombardments. In Fig. [2] the beam is shown striking the 
target mounted on one end of a water-cooled rotating 
arm. A Cu sleeve extended from the cold trap to within 
~1 mm of the target to shield it from condensible va- 
pors. The arm was mounted on a shaft connected via 
a ferrofluidic seal to a computer-controlled servo motor. 
This shaft also carried the water flow for the arm cool- 
ing. A 4 cm by 10 cm copper plate was mounted on the 
opposite end of the arm from the target. This plate con- 
tained 4 apertures with nominal diameters of 1, 2, 3 and 4 
mm, spaced 1.8 cm between centers. The apertures were 
sized by pressing precision-machined steel balls through 
slightly undersized, machined holes, and the aperture ar- 
eas were measured to a precision of ± 0.2% or better 
using a high-magnification traveling microscope. When 
the target was in the bombardment position, the aperture 
plate was directly in front of a Si surface-barrier detector 
as illustrated. Rotating the arm by +180° from its target 
bombardment position placed the 3 mm aperture in the 
beam and put the target in front of the Si a-detector, 
where the 8 B /3-delayed a's were counted. In this con- 
figuration, the beam that passed through the aperture 
was collected in a Faraday cup, shown in Fig. (2J together 
with its electron suppressor ring. A subsequent —180° 
rotation placed the target back in the bombardment po- 
sition. 

A 148 Gd a-source with a moveable shutter was located 
near the Si detector and used to monitor the detector 
gain when the arm was oriented vertically. Several thin 




FIG. 2: Top view of target chamber. 1: shaded beam profile; 
2: LN2 cold trap; 3: target; 4: electron suppressor; 5: Faraday 
cup; 6: ferrofluidic seal; 7: aperture plate; 8: a-detector. The 
3 flanges in the bottom plate (from large to small size) were 
for pumping, heavy- metal shield and the a-detector mount. 

Al fixed shields (not shown in Fig. |2J were placed in the 
chamber and a curved shield was mounted on the aper- 
ture end of the arm to prevent scattered proton beam 
from reaching the a-detector. With this shielding in 
place, no scattered beam was detected above the thresh- 
olds employed in the a-spectrum measurements. 

In order to suppress carbon buildup on the target, the 
chamber was evacuated using a system designed to sup- 
press volatiles. A helium cryopump was used for high- 
vacuum pumping, and sorption pumps for roughing. The 
vacuum system was designed in part for radiation safety. 
7 Be targets were transferred in and out of the chamber 
using a portable Pb-shielded transfer device that allowed 
the target to be mounted on the arm with remote han- 
dling rods. Once the target was mounted, the arm could 
be rotated vertically downward and a heavy-metal shield 
directly below the arm could be raised around the target 
to suppress the 478 keV 7-radiation when the target was 
not in use. A collimated Ge detector was mounted on the 
chamber lid and used to make frequent in situ measure- 
ments of the 7 Be target activity. A large, shielded Nal 
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spectrometer, located on one side of the chamber with its 
front Pb collimator ~ 30 cm from the target, was used for 
7 Be(a, 7) n C measurements of the target energy loss pro- 
files and 19 F(p, a7) 16 resonance measurements of the 
accelerator energy calibration. Horizontal and vertical 
magnetic deflection coils for beam rastering were located 
1.1m upstream of the target, and a magnetic quadrupole 
lens was located 2 m upstream. 



C. Uniform beam-flux technique 

In an in-beam experiment, the yield Y is related to the 
cross section a, the target areal density dN/dA and the 
beam current density dl/dA by the expression 



q J dAdA 



(1) 



where the integral is over the target area and q is the 
charge of a beam particle. If the target areal density is 
constant, and the area of the target is larger than the area 
of the beam, this equation reduces to the usual expression 



Y 



adN 
q dA 



-I, 



(2) 



where I is the total beam current. Eq. is commonly 
used as an approximation in cases where the areal density 
of the target is not constant. However, in this case it is 
difficult to know the error involved, since it depends on 
both the target areal density dN/dA and beam density 
dl/dA nonuniformities, which are often not known. 

On the other hand, if the beam current density is con- 
stant and the target area is smaller than the beam area, 
Eq. reduces to 



q dA 



(3) 



where N is the total number of target atoms. This 
method, which we used, offers several advantages over 
the conventional procedure. We determined N precisely 
by measuring the 7 Be activity, and dl/dA by measur- 
ing the beam transmission through different sized aper- 
tures [2(| . Nevertheless, in any practical experiment the 
beam density can never be made completely uniform. To 
understand possible errors in using Eq. J3J we needed 
to know the density nonuniformities of both the beam 
and target, or more precisely the nonuniformity in the 
product (dN/dA) (dl/dA). We achieved this by relative 
7 Li(d,p) 8 Li measurements as a function of raster ampli- 
tude, as described in Sec. Ill El below. 

We made use of the ti/ 2 = 770 ± 3 ms half- life of 
8 B by bombarding the target and then rotating it away 
from the beam to a shielded location where the 8 B decays 
were detected. A single measurement corresponded to a 
large number of such cycles. Modifying Eq. © to include 



this process, and making all relevant factors explicit, we 
obtain 



cr(E cm ) = 



Y a (E p )F a (E p )f3( s B) 
2(/) p N Bc (t)n/ATT 



(4) 



where E cm is discussed below, E p is the bombarding en- 
ergy, Y a (E p ) is the a yield above a threshold energy, 
F a (E p ) is a correction for the fraction of the a-spectrum 
that lies below the threshold, 4> p is the integrated number 
of beam protons per cm 2 , N-& c (t) is the number of 7 Be 
atoms present at time t of the measurement, and f2 is the 
solid angle of the a-detector. 

The factor /3( 8 B) is the inverse of the timing efficiency 
for counting 8 B decay and is given by 



/3( 8 B) 



Aii[l-e- A ( tl+t2 + t3+t 4)] 
(1 - e At i)[e~ At 2 - e -A(t 2 +t 3 )j ' 



(5) 



where A = 0.693/ii /2 ( 8 B), h = 1.50021 ± 0.00023 s is 
the bombardment period, t 2 = 0.24003 ± 0.00004 s is 
the transfer time to the counting position, t% = t\ is 
the counting period and t 4 = 0.26004 ± 0.00004 s is the 
transfer time to the bombardment position. The arm ro- 
tation was controlled by a hardwired electronic sequencer 
box, and the time periods were measured with a precision 
pulser. Additional tests assured us that the actual arm 
movement and beam on/off periods were consistent with 
these times, with uncertainties similar to or less than 
those listed above, so that /3( 8 B) = 2.923 ± 0.006. 



D. Beam rastering and current integration 

A highly uniform beam flux on target was achieved by 
rastering the beam in x and y using the magnetic deflec- 
tion coils. First, the beam was tuned through the 1 mm 
aperture, with typical transmission 60% or better. Then 
the beam was uniformly rastered over an approximately 
square area, typically about 7 mm x 7 mm, by driving 
the rastering coils with a triangular voltage waveform. 
Incommensurate rastering drive frequencies of 19.03 Hz 
and 43.00 Hz were selected to minimize beating irregular- 
ities in the beam transmitted through a small aperture. 
Pickup coils monitored the time dependence of the mag- 
netic rastering fields. The rastering drive voltages used 
for each cross section measurement were determined from 
aperture scans as described below. 

The arm was electrically isolated from the chamber, 
and biased to +300 V to minimize secondary electron 
losses when the beam struck either the target or the 
aperture plate. The Faraday cup electron suppressor was 
biased to —300 V. The beam currents striking both the 
Faraday cup and the arm were integrated separately, and 
were recorded during both the bombardment and count- 
ing phases. Our notation is: Qt - integrated beam strik- 
ing the target during the bombardment phase; Qa - in- 
tegrated beam striking the aperture during the counting 
phase; and Qc - integrated beam passing through the 
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aperture and striking the cup during the counting phase. 
Just before (after) arm rotation, the beam was rapidly 
swept away from (back onto) the target using one of the 
magnetic deflection coils. 

Our primary integrated beam-charge reference was 
Qc, accumulated during the counting phases of a run 
(many complete arm-rotation cycles). Thus (l/q)Qc/A, 
where A is the area of the 3 mm aperture, is a measure of 
the flux factor cj) p in Eq. (QJ (see below). There are two 
assumptions here: first, that the integrated beam flux 
striking the slightly larger diameter target was equiva- 
lent to Qc/A. Secondly, that the integrated beam flux 
striking the target during the bombarding phases of a 
run was equivalent to Qc/A - that is, as a result of time- 
averaging, the integrated beam flux passing through the 
3 mm aperture and striking the cup during the counting 
period was the same as the integrated beam flux striking 
the target during the bombardment period. The first as- 
sumption, related to beam (and target) uniformity, was 
tested as described in Sec. Ill El below. 

We tested the second assumption by computing the 
factor Q T ■ Qc/(Qc + Qa)- Here Qc/ {Qc + Qa) is the 
fraction of the beam that passes through the 3 mm aper- 
ture, Qt is the total beam striking the target, and thus 
Qt-Qc/(Qc + Qa) measures the (integrated) fraction of 
the beam striking the central 3 mm region of the target. 

In a series of test runs, we found the ratio of the "good 
geometry" (Faraday cup) and "poor geometry" (biased 
arm) beam-flux factors Qc and Qt • Qc /(Qc + Qa) dif- 
fered from unity by 0.02% ± 0.8%, where the uncertainty 
was determined from the run-to-run fluctuations. We 
also found that for short, high-yield 7 Li(d,p) 8 Li diagnos- 
tic runs, the reaction yields showed smaller run-to-run 
fluctuations when normalized to Qt-Qc/(Qc + Qa) than 
when normalized to Qc- Hence we adopted the normal- 
ization (f> p — (l/q)Qr ■ (Qc/A)/(Qc + Qa) based on the 
agreement between the good-geometry and poor geome- 
try results. We made two additional checks on the accu- 
racy of our beam-current integration. We varied the cup 
suppressor bias in the range -300V ± 45V, and found 
<0.5% change in beam current. We also set a limit of 
10~ 4 on the neutral H content of the beam by first tun- 
ing a proton beam onto a LiF target and measuring the 
19 F(p,a7) 16 reaction yield at the 340 and 484 keV reso- 
nances, then sweeping the H + away from the target with 
one of the magnetic deflectors and repeating the measure- 
ment. In this manner we assigned an overall systematic 
uncertainty of ± 0.9% on the integrated beam flux. 



E. Beam and target uniformities 

Figure |3| shows measurements made during the BE1 
experiment to determine the beam and beam-target uni- 
formity. The beam uniformity was determined by mea- 
suring the transmissions through 2, 3 and 4 mm aper- 
tures as functions of (approximately equal) amplitudes 
of the x and y triangular raster waveforms. The top 
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FIG. 3: Top panel: 770 keV deuteron beam transmission ra- 
tios through different apertures, vs. raster amplitude (relative 
units). Bottom panel: 7 Li(d,p) 8 Li yield at 770 keV, normal- 
ized to the integrated beam flux through the 3 mm aperture, 
vs. raster amplitude, measured with the same tune as the 
aperture ratio data. The curves are described in the text. 



panel of Fig. |3| shows measurements made with a 770 
keV deuteron beam, and curves calculated by folding a 
Gaussian beam spot with a rectangular raster distribu- 
tion. The uniformity of the product of the beam and 
target densities was determined by the raster-amplitude 
dependence of the 7 Li(d,p) 8 Li yield from the 7 Be tar- 
get at Ed~ 770 keV, shown in the bottom panel of 
Fig. |21 The curve is a convolution of the target den- 
sity estimated from 7-activity scans, and beam profile 
determined by the transmission ratios, including a fitted 
aperture-target misalignment of 0.5 mm. The behav- 
ior at small raster amplitude was due to this misalign- 
ment. The point at which this yield flattened out deter- 
mined the minimum safe raster amplitude, and is similar 
to the point at which the aperture ratio data flattened 
out. We chose 0.42 as the safe raster amplitude for 770 
keV deuterons, and assigned a conservative ±1% nonuni- 
formity uncertainty here. Aperture transmission curves 
were measured at most proton energies and used, with 
reference to Fig. [31 to determine the raster amplitude for 
each energy and tune so that the beam-target nonunifor- 
mity was < 1%. Independent estimates of the safe raster 
amplitudes were made by folding the target-density dis- 
tribution |25j with beam-flux distributions determined 
from the proton aperture-transmission data, with results 
consistent with the above procedure. The same proce- 
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FIG. 4: Top panel: in situ Ge detector efficiency calibration. 
The curve is a third order polynomial fit. The increase of 
efficiency with energy is due to the absorber. Bottom panel: 
residuals of the fit. 



dure was followed in the BE2 and BE3 experiments. 



F. Be activity measurements 

The 7 Be activity was measured in situ in a "close" 
geometry (target arm vertical and a source-to-detector 
distance of rj 27 cm) by counting the 478-keV 7-rays 
with a 50% efficient Ge detector mounted on the lid of the 
target chamber. An absorber consisting of 8 cm of Al plus 
5 cm of steel reduced the counting rate in the Ge detector 
to 1.5 kHz or less. A cylindrical Pb collimator 6 cm long, 
with a 2 cm hole in the center ensured that the detector 
could not view the small amount of 7 Be sputtered from 
the target and deposited mainly on the cold trap. We 
assumed the accepted 7 Be decay values of t 1 / 2 = 53.12 
± 0.07 d, and BR = 10.52 ± 0.06% for the decay branch 
to the 478 keV level j^. The Ge detector efficiency, 
£478 j was determined using radioactive sources mounted 
on the arm in the same position as the 7 Be target. We fit 
13 lines in the range 276 to 835 keV from 125 Sb, 134 Cs, 
133 Ba, 137 Cs and 54 Mn sources calibrated typically to 
±0.8%(lcr) |H, and obtained ^jv = 1.2. FigureHshows 
the calibration. 

We made an independent check of the Ge detector cal- 
ibration using a second 137 Cs source calibrated indepen- 
dently to ±0.4%(lcr) [2jJ. The result confirmed the cor- 
rectness of our €478 determination at the level of (0.4 ± 
0.8)% . 

In addition, we made a separate determination of the 
BE3 target activity approximately 2 months after the 



end of the BE3 cross section measurements, using the 
Ge detector with a Pb collimator with a 1.6 cm diam- 
eter aperture, no absorbers, a source distance of about 
200 cm and the same calibration sources used in the ear- 
lier measurements. BE3 target activities inferred using 
this calibration and the 53.12 d half-life were 1.5 ± 1.5% 
higher than those determined with the in situ calibra- 
tion. Since these results agree within errors, we used the 
average of these 2 calibrations. 

The number of 7 Be atoms in Eqn.^is given by NB e {t) 
= (3.7 x 10 10 ) A ti/2/ln 2, where the target activity A 
— iV 7 /(3.7 x 10 10 e 4 7 8 i3i?) is given in Curies and iV 7 is 
the number of photopeak 478 keV counts/s in the Ge 
detector. 

There is some evidence that the 7 Be half-life depends 
weakly on host material (see e.g. [!(], El ) ■ We have 
taken the difference between the half-life in Au ( ^53.31 
d) and LiF (53.12 d) as representative of the uncertainty 
due to host material (see Table 1 of j3(j), which intro- 
duces an additional ±0.4% uncertainty in the conversion 
of the measured activities to NBe{t). The total system- 
atic uncertainty in Ns e (t) is ± 1.1%. 

We made a similar, second determination of the BE1 
target activity, in this case 20 months after the BE1 in 
situ measurements. Assuming a 53.12 d half- life, these 
two BE1 activity determinations are in good agreement, 
differing by 0.3 ± 1.9%. Using a 53.31 d half-life the 
difference was 2.7%, suggesting that the 7 Be half- life in 
our host material is close to the accepted value. 

The 7 Be activity measured during the BE3 experiment 
is shown in Fig. [5J As can be seen from the lower panel, 
approximately 15 mCi was lost due to sputtering during 
the cross-section measurements. Compared to a mean 
activity of ^230 mCi during the cross-section measure- 
ments, this is a 7% loss. Hence this activity monitor- 
ing was very important in both the relative and absolute 
cross section determinations. 



G. 'Be target energy-loss profiles 

In previous 7 Be(p,7) 8 B experiments, the energy-loss 
widths of the targets were estimated from the broaden- 
ing of the 41 keV wide 7 Be(p,7) 8 B resonance at E p = 
720 keV [H |H 13, the 12 keV wide E p = 441 keV 
7 Li(p,7) 8 Be resonance or from a calculation using es- 
timated amounts of contaminants 01 ■ 

We determined the complete energy-loss profiles di- 
rectly from measurements of the yield of the narrow 
(r << 1 keV) 7 Be(a,7) 11 C resonance [33J, obtained us- 
ing the large Nal spectrometer. Since the resonance may 
be approximated as a (^-function, and the experimental 
energy resolution was very good, ~ 1 keV, the measured 
7 Be(a,7) 11 C yield directly reflects the energy loss of the 
beam in the target, and may be converted into the corre- 
sponding energy loss distribution for incident protons of a 
given energy using ratios of known proton and a-particle 
energy-loss functions. 
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FIG. 5: Top panel: BE3 target activity versus time since the 
beginning of the (p, 7) measurements. Bottom panel: same as 
top panel but with the decay factor e _At divided out, where 
A = 0.693/53.12 d. Note that the divided activity approxi- 
mately 1700 hr after the end of the experiment agrees with 
the values found immediately after the end of the experiment. 
This demonstrates that the downward slope during the ex- 
periment was due to a loss of target material and not to an 
incorrect half-life. 



The measured (a, 7) profiles were corrected for small 
backgrounds from cosmic rays and the 9 Be(a,n) 12 C re- 
action. In the BE1 experiment, the beam-related back- 
ground was estimated from the observed yield of 4.4 MeV 
7-rays; in the BE3 experiment, by scaling measurements 
of the 7-ray spectrum from 9 Be + a obtained with a 9 Be 
target. 

In the BE1 experiment, the resonance profile was mea- 
sured in the middle and at the end of the (p,7) mea- 
surements. One of these profiles is shown in Fig. 2 of 
ref. |l4j | . The two measured profiles agreed within exper- 
imental error, and the apparent resonance energy repro- 
duced within 1 keV. Thus these measurements show that 
the target profile did not change and there was negligi- 
ble carbon buildup during the (p,7) measurements. We 
found the energy of this resonance to be E a = 1378 ± 3 
keV, in good agreement with the previously determined 
value of 1376 ± 3 keV 1331. The measured BE2 resonance 
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FIG. 6: 7 Be(a,7) u C resonance profiles measured before (top 
panel) , in the middle (center panel) and after (bottom panel) 
the (p, 7) measurements with the BE3 target. The graphs 
show the capture 7-ray yield per /iC of a + beam, per mCi 
of target activity, versus E Q , and the curves are to guide 
the eye. A small nonresonant (neutron) background from 
9 Be(a,n r y) 12 C and a small beam-off background have been 
subtracted from each profile. 



profile was very similar to the El profiles discussed above. 

BE3 target energy-loss profiles were measured three 
times: at the beginning, in the middle and at the end of 
the (p,7) measurements, and are shown in Fig. El Profiles 
2 and 3 are very similar, and differ from profile 1. An 
accident occured during the (p,7) measurements at E p = 
176.5 keV (E cm ~ 139.8 keV) in which the target was ex- 
posed to beam with a small raster amplitude. The close 
similarity of the two profiles measured in the BE1 exper- 
iment, and profiles 2 and 3 in the BE3 experiment shows 
that under normal running conditions both of these tar- 
gets were stable. This provides a strong basis for asso- 
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TABLE I: Target properties. 



Target 


BEl 


BE2 


BE3 


Initial activity (mCi) 


106 


112 


340 


Mean AE a (keV) 


27 


28 


54° 


7 Be:Mo stochiometry 


42:58 


_b 


63:37 



"Profile 2. 

6 No backscattering measurement. 



dating the profile change with the raster accident. The 
(p,7) data taken during the accident were discarded, the 
data measured before the accident were associated with 
profile 1, and later data were associated with profiles 2 
and 3. Proton energy losses inferred from profiles 2 and 
3 were the same within errors. The leading edges of pro- 
files 2 and 3 reproduced to better than AE a = 1 keV, 
and were located at 1 keV higher apparent alpha energy 
than the leading edge of profile 1. This difference was 
within experimental error, indicating that C buildup was 
at most equivalent to AE a = 1 keV. 

The observed profile widths, together with the mea- 
sured 7 Be activity, indicate the presence of target con- 
taminants. From the target-fabrication process |24l l25l ] 
we expect contaminants with masses up to Mo. Analysis 
of 8 B backscattering measurements (see Sec. Ill K|) shows 
that the contaminants are mostly Mo or a material with a 
similar Z. From the profile widths and the known amount 
of 7 Be present, we infer a 7 Be:Mo stochiometry of 42:58 
for BEl and 63:37 for BE3 assuming no other contam- 
inants (see Sec. Ill Kll . These targets were much purer 
than, for example, those of ref. 0. Target properties are 
summarized in Table [I] 

The measured (a, 7) profiles were converted into pro- 
ton energy-loss distributions P{E p ) using the target com- 
position discussed above, and ratios of energy-loss func- 
tions |35l | dE a /dx and dE p /dx, ignoring straggling dif- 
ferences between protons and alphas. Straggling, which 
is significant only in the region of the high-energy tail of 
the energy loss profile, was negligible in our case since 
the high energy tails of our measured profiles decreased 
about a factor of 10 more slowly than expected based 
on straggling calculations |34| , due presumably to target 
nonuniformity. 

Energy-averaged quantities were computed according 

to 

( f(F y _ / f(E P )P(E p )dE p 

{f{Ep)) - JP(E p )dE p ■ (6) 

where f(E p ) is either E p or a(E p ). Energy-averaging 
calculations are discussed further in Sec. IIIII and HVI 



H. Alpha-detector solid angles 

In the BE3 experiment, a-particles from 8 B decay 
were detected in an uncollimated, close geometry with 



a "small" «139 mm 2 20/im Si surface-barrier detector, 
that we label "S" , and, in separate measurements, with 
a "large" «416mm 2 33/zm detector, that we label "L". 
We inferred the solid angle ft of each of these detectors 
from the counting rate of a custom-made 148 Gd a-source 
deposited on a Mo target backing of the same design as 
was used for the 7 Be target. The source activity was cal- 
ibrated with a ~450 mm 2 35/x detector collimated to an 
area of 246.0 ± 0.3 mm 2 , and a source-to-collimator dis- 
tance of 57.00 ± 0.10 mm, corresponding to a geometrical 
solid angle ft ca i = 0.0744 ± 0.0003 sr. The efficiency of 
the calibration detector for events lying within the colli- 
mator acceptance was checked by measurements with dif- 
ferent size collimators and different detectors. From the 
measured a-source counting rates we determined ft/ft ca i 
with a statistical precision of ± 0.2%. We applied (0.6 
± 0.6)% and (0.3 ± 0.3)% corrections to the deduced 
solid angles for the S and L detectors, respectively, to 
account for a measured (0.03 ± 0.03) mm difference in 
the distance from the detector to the BE3 target and to 
the 148 Gd source. The results are ft/ft ca i = 24.17 ± 0.39 
and 18.78 ± 0.13, and hence ft = 1.798 ± 0.030 sr and 
1.397 ± 0.011 sr for the S and L detectors, respectively. 

To understand better the accuracy of these solid- 
angle measurements, we measured the L-detector effec- 
tive solid angle as a function of distance using a pre- 
cision translation stage. Several different sources, in- 
cluding commercially- produced 148 Gd and 241 Am sources 
were employed. Based on the observed spread of these 
measurements, we conservatively assigned an additional 
scale- factor (common mode) error of ± 1.5% to our de- 
terminations of ft for both the S and L detectors (see 
Table CU). 

This method had several improvements over the BEl 
experiment 0|. 1) It eliminated the 7 Li(d,p) 8 Li solid- 
angle calibration which depended on calculated correc- 
tions for the portion of the 8 Li a-spectrum lying below 
the experimental threshold. 2) The smaller a-detector 
solid angles helped reduce the solid-angle uncertainty. 3) 
The thinner detector minimized corrections for the por- 
tion of the 8 B a-spectra lying below the experimental 
threshold. 



I. Beam-energy calibration 

The accelerator-energy calibration in the BEl exper- 
iment was determined by measuring the 19 F(p,a7) O 
resonances at 340.46 ± 0.04, 483.91 ±0.10 and 872.11 ± 
0.20 keV [3(| using a thick LiF target. The measured 
thick-target yield curves were fitted with the integral of 
a Lorentzian folded with a Gaussian beam energy resolu- 
tion. Multiple measurements of each resonance, in which 
the beam tune and steering were varied, were used to 
estimate the systematic error. 

In the BE3 experiment, we remeasured these same res- 
onances, as well as the 19 F(p,a7) 16 resonance at E p 
= 223.99 ± 0.07 keV , and the 7 Be(a,7) n C resonance 
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FIG. 7: Accelerator energy calibration constant k vs. proton 
energy. Points: measured values; solid line: spline fit. The 
error bars were determined from their scatter relative to the 
fit. The 7 Be(a,7) 11 C resonance measurement (at an equiva- 
lent proton energy of 5471 keV) indicated that k was constant 
above E B = 1000 keV. 



which we take to be located at 1377 ± 2 keV from the 
mean of our resonance energy determination in the BE1 
experiment and the determination given in ref. [33j. A 
spline curve was fitted through nine measurements of 
these five resonances and used to determine the accel- 
erator calibration constant fc, shown in Fig. Q where 



EL 



kf 2 



1 + E p /2m p c 2 



(7) 



and / is the magnet NMR frequency. The BE2 and BE3 
measurements were made under similar accelerator con- 
ditions, and hence we used the same accelerator calibra- 
tion for both of these experiments. 



J. Alpha spectra from 7 Be(p,7) 8 B 




E a (keV) 



FIG. 8: Alpha spectra measured with the BE3 target and the 
S and L detectors, at different proton energies, as indicated. 
The curves are fitted TRIM calculations. 



Fig. [S] shows representative a-spectra measured with 
the BE3 target. The energy scale was determined using 
a 148 Gd source to monitor the detector gain and a preci- 
sion pulser to monitor the zero. Both gain and zero were 
checked frequently during the experiment, and found to 
be stable, the gain to ± 0.2% and the zero to ± 3 keV. At 
low energies the a-spectrum was obscured by background 
from pileup of Compton electrons produced by the in- 
tense 478 keV radiation from the target, and at even 
lower energies by electronic noise. The a-spectra were 
summed above a software threshold of 650 keV (900 - 
930 keV) in the S (L) detector, and corrected for beam-off 
background of 0.6 ± 0.1 counts/hr (averaged between the 
S and L detectors), due primarily to a-radioactivity in 
the chamber materials. A stainless-steel shield between 
the chamber wall and the Si detector reduced this back- 
ground by a factor of 2 relative to the (small) level ob- 



served in the BE1 experiment. This background correc- 
tion was 5%, 2% and 1% at the 3 lowest proton energies, 
and less than 0.3% at higher E p . Beam-on background 
was checked at several energies and found to be negligi- 
ble. Separate measurements with a LiF target ensured 
that any deuteron beam contamination (that could con- 
taminate the /3-delayed a-spectrum by the 7 Li(d,p) 8 Li 
reaction) was negligible. 

To compute the total yield, the data integrated as 
described above must be corrected for the fraction 
FaiEp) -1 of the spectrum that lies below the detector 
threshold. This was done using the TRIM Monte Carlo 
(MC) code [35| to model the implantation depth of the 
8 B ions and the energy loss of the emitted ev's. The 
a-energy spectrum was taken from the thin-target spec- 
trum of ref. |37| . Calculations were performed in which 
a) the opening angle of the emitted a's was restricted to 
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the geometrical acceptance of the detector, and b) the a's 
were randomly emitted into 2tt. In case b), which allows 
for multiple scattering into and out of the detector ac- 
ceptance, additional events were found at low E a < 200 
keV due to large-angle multiple scattering. However, the 
F a (Ep) values predicted by the 2 methods agreed within 
the MC precisions of ± 0.2%. The curves shown in Fig. El 
were computed with method b), for which 



F a (Ep) 



QN a (E a )dE a /N tot ' 



(8) 



where f2 is the geometrical solid angle, N a (E a )dE a is 
the integral of all MC events above the detector thresh- 
old, and Ntot is the total number of MC events. 

The S detector spectra (except for the highest energy 
point at E cnl — 1754 keV) were fitted with a fixed en- 
ergy calibration, varying only the counting-rate normal- 
ization. For the 1754 keV data point and for most of 
the L-detector data, we had to vary the zero offset to 
fit the low-energy side of the spectra. The resulting fit- 
ted offsets (typically 20 - 30 keV with a maximum of 88 
keV) are larger than can be accounted for by the mea- 
sured offsets. They are also larger than the 1-2 keV 
expected, on average, from (3 - a summing, and suggest 
a failure of the TRIM calculation. A similar problem 
occured in comparisons of TRIM calculations with mea- 
sured 7 Li(d,p) 8 Li spectra. As a result of these difficulties, 
and also problems in the BE2 experiment understanding 
the ratio of 7 Li(d,p) 8 Li yields to a-source count-rates for 
different detector solid angles, we assigned a conservative 
uncertainty of ± 30% on the correction 1 — F a (E p ). 

For the S (L) detector, F a (E p ) ranged from 1.004 to 
1.006 (1.023 to 1.037 ) below the resonance, and 1.012 to 
1.019 (1.054) above the resonance . 

Threshold correction factors for the BE2 experiment 
are discussed in Sec. IIVI below. 



K. Backscattering measurements 

Weissman et al. [38| pointed out that previous 
7 Be(p,7) 8 B experiments suffered from unknown losses of 
8 B due to backscattering out of the target, and 8 Li losses 
when 7 Li(d,p) 8 Li was used for absolute cross section nor- 
malization. Substantial backscattering losses may occur 
when a high-Z target backing is used or if there are high- 
Z contaminants in the target. Proton backscattering fol- 
lowed by 8 B production is not important. 

We have made the only measurements of 8 B backscat- 
tering in the 7 Be(p,7) 8 B reaction, by modifying our ap- 
paratus as shown in Fig. |U| We installed the 7 Be target 
in a fixed mount in place of the Faraday cup, mounted 
large-area Cu catcher plates on both ends of the ro- 
tating arm, and installed the L-detector on the down- 
stream side of the arm. During target bombardment, 
the proton beam passed through a 4 mm aperture in 
the center of the catcher plate before striking the tar- 
get, and backscattered 8 B's were caught on the catcher 




FIG. 9: Backscattering apparatus. 1 and 2: catcher plates; 3: 
fixed target and water-cooled mount; 4: a-detector. 



plate. The arm-rotation time sequence was the same as 
in the Be(p,7) 8 B measurement. The arm and the tar- 
get were both biased to +300V. Because of secondary 
electron crosstalk between the target and the arm, only 
the total beam current could be measured reliably in this 
setup. Hence beam tuning and beam transmission mea- 
surements through the 4 mm aperture were done with 
the Faraday cup in place, after which the 7 Be target was 
installed and the backscattering measurements were car- 
ried out. The fraction of the beam passing through the 
aperture was checked at the end of each measurement 
by re-installing the Faraday cup. Backscattering mea- 
surements were made at E cm — 626 and 1200 keV. The 
efficiency for catching the 8 B on the plates and count- 
ing the subsequent /3-delayed a-particles was computed 
using TRIM to estimate the backscattered 8 B angular 
distribution. Figure EH shows the measured backscat- 
tering probabilities for the BE1 and BE3 targets, along 
with TRIM calculations for several assumed target com- 
positions. We assumed uniform targets composed of 7 Be, 
C and Mo, where C is representative of low-Z contami- 
nants and Mo is representative of high-Z contaminants 
(see Sec. Ill C|l . The amount of 7 Be in the target is fixed 
by the target activity, and the amount of contaminants 
is fixed by the (a, 7) resonance profile width in excess of 
the width expected for pure 7 Be (see Sec. lIICll . 

The qualitative shapes of these curves are easy to un- 
derstand. For pure C contaminant, the backscattering 
is primarily from the higher-Z Mo backing. Although 
the (Rutherford) backscattering cross section rises as the 
energy drops, at low bombarding energy the backscat- 
tered 8 B ions have insufficient energy to escape from the 
target layer. For pure Mo contaminant, the backscat- 
tering probability is higher and extends to lower ener- 
gies, since the backscattering may now occur in the tar- 
get layer. The backscattering probability from BE3 is 
smaller than from BE1 because this target was thicker 
which suppressed backscattering from the Mo backing at 
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FIG. 10: Measured backscatting probabilites (errors include 
statistics and systematics) and TRIM calculations. The 
7 Be:C:Mo atom-number ratios assumed in the TRIM calcu- 
lations are, in descending order: top panel - 63:0:37, 58:8:34, 
57:13:30; bottom panel - 42:0:58, 38:19:43, 36:31:33 and 
19:81:0. 
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FIG. 11: BE3 S (circles) and L (squares) S-factor data mea- 
sured below the resonance. Curve: best DB fit. Error bars 
include statistical and varying systematic errors. 



III. BE3 RESULTS 
A. Data 

Figure ITU shows the data from the BE3 experiment. 
The cross sections were converted into S-factors using 
the relation 



S 17 (E cm ) — &(E cm )E c 



(9) 



where the Gamow energy Eg = (27raZiZ2) 2 /ic 2 /2 = 
13799.3 keV, and E cm values were computed as described 
below. This procedure corrects for energy averaging, so 
that fitting experimental Sn(E cnl ) values without explic- 
itly including energy averaging is equivalent to fitting 
the measured cross sections including energy averaging. 
TablelTTllists the uncertainties in the experimental Si7(E) 
values. 



low energies. 

As shown in Fig.^J| the best fit curves are 7 Be:C:Mo 
= 38:19:43 for BE1, and 7 Be:C:Mo = 58:8:34 for BE3. 
These contaminant compositions are similar, which is 
expected since the target fabrication process was the 
same. We note that a pure Mo contaminant is con- 
sistent with both target compositions, within errors. 
The backscattering probability can be affected by target 
nonuniformity. Even though BE1 was significantly less 
uniform than BE3, as can be seen from the 7 Be(a, 7) n C 
profiles, the agreement between measured and calcu- 
lated backscattering probabilities with similar contam- 
inant compositions suggests that nonuniformity did not 
play an important role. 

Based on Fig. ^3 we made constant backscattering 
corrections to our measured cross sections of 0.4 ± 0.1% 
for BE3 and 1.0 ± 0.5% for BEL 



B. E cm values 

We computed the effective energy E cm of the data 
points as follows. First, we computed (<r(E cnl )) using the 
appropriate 7 Be(a,7) 11 C profile (see Sec. Ill where 
o'(Ecm) is given by Eq. © with Si7(E cm ) set equal to 
a constant. Here we assumed a pure Mo contaminant 
which is consistent with the backscattering data - see Sec. 
Ill Kl Then we solved the equation (<r(E cm )) — <r(E cnl ) 
for E cla . Our neglect of the variation of Sn(E cm ) over 
the target thickness is a good approximation everywhere 
except near the resonance. Our analysis of the BE1 data, 
which included the resonance, took into account the ef- 
fect of the resonance on the calculated -E C m's - see Sec. 
IfVI Away from the resonance, E cm and (E cm ) values are 
similar everywhere except at the lowest energies; for ex- 
ample, at E p = 149.9 keV, the energies E cm = 113.9 keV 
and (-E^m) = 115.6 keV correspond to a 6% difference in 
Sn(E cm ). 
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TABLE II: Percent uncertainties in S17 from the renormalized 
BE3 (S and L) data. 



Statistical errors 


1.3- 


-4.0 


Varying systematic errors: 






proton energy calibration 


0.0 - 


-0.7 


target thickness 


0.0- 


- 1.6 


target composition 


0.0- 


-0.7 


a-spectrum cutoff 


0.1 - 


- 0.7(5) 




1.0 - 


- 1.8(L) 


kJVjdlC Id/L bUl CI L Ul O • 






beam-target inhomogeneity 


1.0 




integrated beam flux 


0.9 




7 Be target atom number 


1.1 




solid angle 


1.5 




backscattering 


0.1 




timing cycle 


0.2 




Total scale factor error 


2.3 






50 100 150 200 250 300 350 400 



FIG. 12: BE3 S and L S-factor data renormalized to the best- 
fit Si 7 (0) = 22.1 eV b, otherwise as in Fig. [HI 



IV. BE1 AND BE2 RESULTS 



C. Determination of Si7(0) 

Because the experimental data must be extrapolated 
to low energy to determine the astrophysical S-factor, 
it is best to fit data as low in energy as possible, com- 
mensurate with good experimental precision. Above the 
Ml resonance at E cm = 630 keV, different 7 Be(p,7) 8 B 
cross-section calculations deviate substantially from one 
another. Below the resonance, the 7 Be(p,7) 8 B reaction 
is predominantly direct capture, and becomes increas- 
ingly extranuclear and hence less model-dependent with 
decreasing E p . 

The 7 Be(p,7) 8 B calculation that fits experimental data 
best over a wide range is the cluster-model theory of Dc- 
scouvemont and Baye 39] (DB) (see Sec. IV Al below) . 
FigurelTTlshows our fit of the scaled DB theory to all our 
data with E cm = 116 to 362 keV. The DB theory shown 
here does not contain contributions from either the 1 + 
resonance near 630 keV or the 3 + resonance near 2200 
keV. For these low-energy data the contribution of the 
1 + resonance is less than 0.4% and the contribution of 
the 3 + resonance is completely negligible based on our 
resonance fits discussed below. The resulting Si7(0) val- 
ues are 22.35 ± 0.41 eV b and 21.95 ± 0.34 eV b for the 
S and L detector data, respectively. These values agree 
well within the quoted errors, which include only statisti- 
cal and relative (independent) systematic contributions. 
Our best value for Si7(0) is the weighted average of these 
2 results. Including the common-mode scale-factor error 
of ±2.3% from Table UH we obtain 

Si 7 (0) = 22.1 ± 0.6(expt) eV b , (10) 

where the error includes all contributions other than 
the theoretical extrapolation uncertainty discussed be- 
low. The curve shown in Fig. ^2 corres P on ds to the DB 
calculation scaled to our best fit value for Si7(0) as given 
in Eq. I|10[l above. Fig. 1121 shows these same data renor- 
malized to the best fit value for Si7(0). 



Results for our BE1 experiment were reported in 
ref. 0. Figure ED an d Table show these results, 
renormalized to our best value for Si7(0) as given in 
Eq. (|10fl . The error bars have been increased relative to 
those shown in Fig. 3 of ref. 0], because of the increased 
F a (E p ) uncertainty discussed above; otherwise, the data 
are the same as in ref. |14| . Note that our previously pub- 
lished value, 517(0) = 22.3 ± 0.7(expt) eV b 14], is in 
excellent agreement with our new value presented above. 

In the BE2 experiment, data were taken in a close ge- 
ometry with the L detector, over the range E cm = 876 
to 2459 keV. These data are shown in Fig. ^] and ITSI 
where they are plotted, along with our BE1 and BE3 
data with a common absolute normalization determined 
by our BE3 results (Eq. The BE2 data have rela- 

tively large systematic errors due to the large threshold 
correction factors F a (E p ) that range from 1.091 to 1.287, 
and the corresponding ± 30% uncertainty on 1 — F a (E p ). 



A. 1 + resonance 

In the BE1 experiment, data were taken from E cnl = 
186 to 1200 keV, including detailed measurements over 
the 630 keV Ml resonance. E cm values were determined 
by fitting the cross-section data, including averaging over 
the target profile, with the sum of the nonresonant DB 
cross section plus an incoherent Breit-Wigner resonance 
(see below), and then solving the equation a = cr(E cul ) 
for E cm , where a is the measured cross section and 
c(-Bcm) is the fitted unaveraged cross section at a given 
proton bombarding energy. This procedure removes the 
effects of energy averaging, as discussed above. Measured 
cross sections were converted to S- factors using Eq. (0. 
We checked that fitting our Si7(i? cm ) data without en- 
ergy averaging, gave the same parameters as the original 
fit to the cross section including energy averaging. 
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FIG. 13: Top panel: BE1 data normalized to Si 7 (0) = 22.1 eV 
b. Solid curve: best-fit DB plus a fitted 1 + resonance; dashed 
curve - DB only (see text). Inset: resonance region. Bottom 
panel: solid squares - BE3 S data; solid circles - BE3 L data. 
The solid curve was calculated with 1 + resonance parameters 
determined from fits to the BE1 data, and the normalization 
was determined by fitting the BE3 data with E cm < 362 keV. 



The cross section was fitted with 



cr(E cm ) = Cia DB {E cm )+ 



Co 



r 'p(-^crn)^" 7 (-^cm) 



E cm (E cm -E a )2 + r p (E cm )y4 

_ (11) 

where C% « 0.7 is a fitted scaling factor, obb (E C m) is 
the DB cross section (with the 1 + resonance removed), 
C 2 = 3nX 2 E cm /8 , T p (E cm ) = T p (E )Pi(E cni )/P 1 (Eo), 
Pi(E cm ) is the £=1 Coulomb penetrability evaluated at 
R = 3.65 fm, T 7 (^ cm ) = T 7 (E )E^/El E. ( = E cm + Q 
and Q = 0.137 MeV. 

Table IIIII shows our center-of-mass 1 + resonance fit 
parameters together with those of refs. 0, ^| and the 
recent elastic scattering results of ref. |40|. Descouve- 
mont and Baye [3^ predict the lowest 1 + resonance at 
E cm ~ 0.2 MeV. Scaling by the experimentally measured 
energy, they calculate F p (E ) w 59 keV and T 1 (Eq) = 33 
meV (assuming pure Ml), in reasonable agreement with 
experiment. 
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FIG. 14: Cross-sections from the BE1, BE2 and BE3 exper- 
iments. The BE2 data are shown as open squares, and the 
BE1 and BE3 data are shown with the same symbols as in 
Fig. 1131 Solid curve: best-fit DB plus fitted 1 + and 3 + reso- 
nances; dashed curve: DB only; lower solid curve: 1 + and 3 + 
resonance contributions; as in Fig |13l 
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2500 



FIG. 15: Measured and calculated S-factors corresponding to 
the cross sections shown in the previous figure. Solid curve: 
best-fit DB plus fitted 1 + and 3 + resonances; dashed curve: 
DB only; lower solid curve: 1 + and 3 + resonance contribu- 
tions. 



U. a resonance 

Our data in Fig. ^| show clear evidence for the low- 
est 3 + resonance at E cm ~ 2200 keV. The DB calculation 
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TABLE III: 1 + resonance parameters 



Parameter 


Present work 


Ref. [9] 


Ref. [15] 


Ref. [40] 


E (keV) 


630 ± 3 


632 ± 10 


633 


634 ± 4 


T P (E ) (keV) 


35.7 ± 0.6 


37 ± 5 


35 ± 3 


31 ± 4 


r 7 (£ ) (meV) 


25.3 ± 1.2 


25 ± 4 


25 ± 2 





TABLE IV: 3 + resonance parameters 



Parameter 


Constrained fit 


Unconstrained fit 


E (keV) 


2183" 


2100 ± 60 


T P (E ) (keV) 


350° 


510 ± 270 


r 7 (£ ) (meV) 


150 ± 30 


180 ± 70 



ref. El . 



shown in Figs. El]-EJ3does not include this 3 + resonance. 
To understand better the 7 Be(p,7) 8 B cross section above 
E cm = 1200 keV, we have fitted our BE2 + BE3 data 
with the DB calculation plus adjustable 1 + and 3 + res- 
onances. The 1 + parameters were fixed to our fit results 
quoted in Table lllll and the 3 + resonance was fitted with 
a formula similar to that given in Eq. (|llf> . in which the 
constant Ci was multiplied by the factor 7/3 to account 
for the J = 3 resonance angular momentum. This for- 
mula neglects f-wave capture and El decay, which should 
be good approximations. We also neglect the 4 He + 3 He 
+ p channel and (here as in the 1 + resonance analysis) 
proton inelastic scattering to the first excited state of 
7 Be. 

The fit results are shown in Table IIVI The uncon- 
strained fit parameters, E = 2100 ± 60 keV, T p (Eq) = 
510 ± 270 keV and r 7 (£ ) = 180 ± 70 meV, are not 
well determined. The resonance energy and width agree 
with the more precise values 2183 ± 30 keV and 350 ± 40 
keV, respectively, compiled in |4l| . Constraining the res- 
onance energy and width to the values from |4l| results in 
the fits shown in Figs. IbH and ITSI for which Tj(E ) = 150 
± 30 meV. Descouvemont and Baye [3!| also calculated 
this 3 + resonance, and found Eo ~ 2800 keV. After ad- 
justing their resonance energy to agree with experiment, 
they obtained T p (E ) = 530 keV and T 7 (£;o) = 45 meV. 

The DB calculation together with fitted 1 + and 3 + res- 
onances provides a reasonable description of the data up 
to E cm = 2500 keV, although the constrained fit indicates 
that the DB "background" underneath the 3 + resonance 
is not quite correct. Including the fitted 3 + resonance 
has very little (< 1%) effect below 1500 keV; hence we 
do not include the fitted 3 + resonance elsewhere in this 
paper. 



TABLE V: S17 data normalized to the best-fit BE3 results, 
and lcr errors from this work. <r a t a t - statistical error, a va , ly - 
varying systematic error, all in eV b. An additional common- 
mode error of 2.3% applies to each point. 





















Sit 


Cstat 


<^~vary 


E crn 


S17 


fstat 




BE3 S 


BE3 L 


184.3 


19.8 


0.4 


0.2 


110. D 


20.6 


0.8 


n a 


219.8 


19.3 


0.4 


0.2 


loy.o 


19.2 


0.6 




255.4 


19.4 


0.3 


0.1 


lOI .u 


18.9 


0.5 


n 9 
u.z 


277.5 


20.1 


0.3 


0.1 


9^ *} 
ZOO.o 


19.5 


0.4 


n 9 

U.Z 


326.4 


20.7 


0.4 


0.1 


977 ^ 


19.9 


0.4 


n 9 
u.z 


361.9 


20.1 


0.3 


0.1 


Q9A A 


19.9 


0.3 


n 9 
u.z 


871.2 


24.3 


0.3 


0.1 


ooi.y 


20.9 


0.4 


n ' -i 


999.5 


24.6 


0.3 


0.1 


S271 A 


25.0 


0.4 


n a 


1099.8 


25.6 


0.3 


0.1 




BE 2 




1200.1 


26.4 


0.6 


0.1 




24.4 


0.2 


U . 


1754.1 


30.8 


0.8 


0.2 


1UU1.D 


24.5 


0.2 


u . 


BE1 


14UO.O 


27.4 


0.4 




185.6 


19.3 


0.6 


0.4 


id / y.4 


29.3 


0.5 


u . 


221.3 


19.3 


0.5 


0.3 


i q^i n 
lyoi.u 


34.8 


0.6 


U.o 


257.0 


19.4 


0.6 


0.3 


91 QA 7 


43.2 


0.8 


n o 
u . y 


293.5 


19.6 


0.5 


0.3 


Z40O.0 


40.6 


0.7 


n o 
u.y 


294.4 


20.2 


0.6 


0.3 




BE1 






328.2 


20.3 


0.4 


0.3 


639.4 


89.3 


1.4 


1.7 


363.8 


20.3 


0.4 


0.3 


649.2 


60.9 


1.2 


1.1 


408.1 


20.6 


0.3 


0.3 


658.7 


47.5 


0.9 


0.9 


461.3 


21.1 


0.4 


0.3 


679.1 


33.2 


0.9 


0.6 


496.7 


22.0 


0.3 


0.4 


699.4 


27.0 


0.7 


0.5 


528.6 


22.8 


0.6 


0.4 


750.7 


25.5 


0.7 


0.5 


558.8 


25.4 


0.7 


0.4 


820.7 


23.9 


0.5 


0.5 


589.0 


34.9 


0.8 


0.6 


876.3 


24.3 


0.2 


0.5 


599.7 


44.8 


1.1 


0.8 


876.3 


24.0 


0.4 


0.5 


609.4 


56.9 


1.4 


1.0 


1002.3 


24.2 


0.2 


0.6 


619.6 


86.8 


1.6 


1.6 


1102.8 


25.2 


0.3 


0.6 


633.3 


103.7 


1.1 


1.9 


1203.2 


25.4 


0.3 


0.7 



V. DISCUSSION 
A. Extrapolation uncertainty 

The DB calculation with an empirically-fitted 1 + reso- 
nance provides the best description of the energy depen- 
dence of the 7 Be(p,7) 8 B cross section in the range E cm < 
1200 keV, although it underestimates the data slightly in 
the region near 1200 keV. 

It seems reasonable that the DB calculation, which 
fits the 7 Be(p,7) 8 B data best over a wide energy range, 
should be most reliable for extrapolating to energies of as- 
trophysical interest. Therefore, we base our central value 
for Si7(0) on the DB extrapolation of our low-energy 
data. However, it is possible that calculations that do 
less well at high energies may be acceptable at low ener- 
gies. We estimated the extrapolation uncertainty by fit- 
ting 12 published calculations ^2] to our BE3 data with 
Ecm < 400 keV. The results, shown in Fig. exhibit a 
total spread of 2 eV b (out of 22 eV b) at zero energy, and 
an rms deviation of ± 0.6 eV b which we adopt as the 
theoretical extrapolation uncertainty. Hence our best-fit 
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TABLE VII: Experimental Si7(0) values and uncertainties in 
eV b determined by our DB fits to published data. 
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FIG. 16: Fits of 12 different theories to the BE3 data below 
the resonance - see 14a. 



TABLE VI: 517(20) and S 17 (0) (in eV b) from fitting our data 
with _E cm < 362 keV with different models, as in Fig. 1161 and 
ref. IH. 



Model 


5-17(20) 


517(0) 


Nunes 


20.8 


21.4 


Johnson 


20.5 


21.2 


Bennaceur 


21.5 


22.2 


Barker B80 


20.7 


21.2 


Barker Bl 


21.8 


22.6 


Barker B2 


21.1 


21.8 


Csoto C2B 


21.7 


22.0 


Csoto C8B 


21.8 


22.1 


Jennings r c = 2.4 fm 


22.0 


22.8 


Jennings r c = 1.0 fm 


21.1 


21.8 


Typel 


20.3 


20.8 


Descouvemont 


21.4 


22.1 



result is 

Si 7 (0) = 22.1 ± 0.6(expt) ± 0.6(theor) eV b. (12) 

Our theoretical error estimate of ± 0.6 eV b is some- 
what larger than the value quoted in ref. |l4j because we 
now include the Typel calculation. It is also considerably 
larger, and hence more conservative than the ± 0.2 eV 
b uncertainty recommended by Jennings Note that 
the theoretical (extrapolation) uncertaintly is as large as 
the experimental uncertainty so that additional theoret- 
ical work to reduce the extrapolation uncertainty would 
be very valuable. 

In solar model calculations, the value of Si7(0) along 
with the derivatives S' 17 (0)/Si 7 (0) and Si' 7 (0)/Si 7 (0) 
are used to compute the 7 Be(p,7) 8 B reaction rate in the 
sun. However, Jennings has pointed out that the 
derivatives vary significantly among the different theo- 
ries, and also differ from the "best" values given in ref. Q . 



Fit Range 


< 425 keV 


< 1200 keV 


Experiment 


Value 


Error 


Value 


Error 


Filippone 


20.7 


2.5 


19.4 


2.2 


Hammache 


20.1 


1.3 


19.4 


1.1 


Hass 






20.4 


1.1 


Strieder 


18.8 


1.8 


18.1 


1.6 


Baby 


20.8 


1.3 


21.2 


0.7 


This work 


22.1 


0.6 


22.3 


0.6 


<425 keV best fit 21.4 
<1200 keV best fit 


0.5 


21.1 


0.4 



He |43j argues that Si 7 (20) should be used in solar model 
calculations, which avoids the need for derivatives since 
20 keV is near the center of the Gamow window. Using 
Si 7 (20) instead of Si 7 (0) also avoids the need to extrapo- 
late theoretical cross section calculations to zero energy. 

For these reasons, we also quote our best-fit result for 
the S-factor at 20 keV: 

S i7 (20) = 21.4±0.6(expt)±0.6(theor) eV b. (13) 



Table IVII displavs our Si 7 (0) and Si 7 (20) values obtained 
from the different theoretical extrapolations. There is 
a 3% variation in the ratio Si 7 (20)/Si 7 (0) among these 
theories, which is surprisingly large. 



B. Comparison with other direct experiments 

We compare the results of all direct experiments by 
fitting the DB theory to published data in two different 
energy ranges: E clrL < 425 keV, and E cm < 1200 keV. We 
made a substantial effort to ensure accuracy in these com- 
parisons by obtaining data from primary sources when- 
ever possible, and by fitting the data ourselves. 

In the low-energy range, the 1 + resonance contribution 
may be neglected (it is « 1% of the direct contribution 
at E cm = 425 keV and drops rapidly with decreasing en- 
ergy) , and the theoretical uncertainty is minimized. The 
experimental uncertainty due to the a-threshold correc- 
tion is also minimized, and the high energy tail of the 1 + 
resonance is avoided. On the other hand, some experi- 
ments do not have good precision at low energies (none 
are as good as the present study), which motivates our 
wider-range comparison. In the wide-range fits we in- 
cluded the 1 + resonance with parameters fixed from the 
fit to our data and excluded data close to the resonance. 
Care was taken to separate common-mode (scale factor) 
errors from other errors. 

We renormalizcd all published data that used the 
7 Li(d,p) 8 Li normalization to cr[ 7 Li(d,p) 8 Lil = 152 ± 6 
mb, the average of the results quoted in 0, HU for the 
cross section at th e p eak of the broad 780 keV reso- 
nance. The results 44] for both fitting ranges are shown 
in Fig. [T7|and Table IVTT1 
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TABLE VIII: Experimental Si7(20) values and uncertainties 
in eV b determined by our DB fits to published data, other- 
wise the same as Table IViTI 

Fit Range < 425 keV < 1200 keV 



Experiment" 


Value 


Error 


Value 


Error 


Filippone 


20.0 


2.4 


18.8 


2.2 


Hammache 


19.4 


1.2 


18.8 


1.0 


Hass 






19.7 


1.0 


Strieder 


18.1 


1.7 


17.5 


1.5 


Baby 


20.1 


1.3 


20.5 


0.7 


This work 


21.4 


0.6 


21.5 


0.6 


<425 keV best fit 
<1200 keV best fit 


20.6 


0.5 


20.3 


0.4 



"Individual Si7(20) values were computed from the original data 
before rounding. 
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FIG. 17: 517(0) values determined from our DB fits to pub- 
lished data from direct experiments. Bottom panel (top 
panel): fits to data with E cul < 425 keV (1200 keV). The hor- 
izontal solid and dashed lines indicate the mean values and 
uncertainties determined from fitting the data of Filippone 
and more recent experiments. 



Results from [U IE 0, ma y suffer additional error 
from 8 B and 8 Li backscattering losses. We estimated 
these losses for Filippone's experiment Ql using TRIM 
and the target compositions described in [45|]. We calcu- 
lated both the loss of 8 B from 7 Be(p,7) 8 B and the loss 



of 8 Li from the 7 Li(d,p) 8 Li reaction, the latter having 
been used by for one of the two absolute cross-section 
determinations. We find smaller corrections by a fac- 
tor of 2 compared to those published by Weissman [38| . 
Our calculated corrections for the weighted average of 
Filippone's normalizations, approximately —2% to —4% 
depending on proton energy, are sufficiently small that 
we ignore them. Hammache et al. |10j applied calcu- 
lated backscattering corrections to their data; Strieder 
et al.0] used a low-Z backing and assumed negligible 
losses; and Baby et al. used an implanted target for 
which losses should have been negligible. We made no 
corrections to any of these published data for backscat- 
tering losses. 

We combined all the direct S-factor data shown in 
Fig. 1171 except for the older pioneering experiments of 
Parker Kavanagh |(j and Vaughn Q. The results 
for E cm < 425 keV are, Si 7 (0) = 21.4 ± 0.5 eV b, ^jv 
= 1.2 {u = 4), and for E cm < 1200 keV, Si 7 (0) = 21.1 
± 0.4 eV b, x 2 /" = 2.1 {v = 5). These best fit values 
and uncertainties are listed also in Table IV11I The fit 
to the low-energy region has a good x 2 , with P(x 2 ,v) 
= 30%. The wide-range fit has P(x 2 ,v) = 0.06, which 
is not unreasonable, though it suggests that some of the 
experimental uncertainties may be underestimated. The 
good agreement between the Si 7(0) values determined 
from the 2 different fit ranges demonstrates that the en- 
ergy dependence of the 7 Be(p,7) 8 B cross section is well- 
understood and similar for all the direct experiments, 
and that the Si 7 (0) value determined here is insensitive 
to the fit range. We show in Table IVTTT1 th e best-fit DB 
results for Si 7 (20), analogous to Table IvTlI for Si 7 (0). 

As argued above, low-energy data provide the most re- 
liable basis for determining Si7(0). Hence our best value 
determined from direct experiments is 

S i7 (0) = 21.4±0.5(expt) ± 0.6(theor) eV b, (14) 

where we have taken the theoretical extrapolation error 
from Sec. IVAl 

The energy dependence of the 7 Be(p,7) 8 B cross sec- 
tion is also well-determined by our measurements, since 
we have minimized all of the important uncertainties 
here. Our data have small statistical uncertainties and 
small point-to-point scatter, both above and below the 
1 + resonance. Our small area target - uniform beam flux 
technique allowed us to determine the number of beam- 
target interactions reliably for different beam energies 
and beam focussing conditions. Our frequent in situ ac- 
tivity measurements determined the number of 7 Be tar- 
get atoms present during each (p,7) measurement, avoid- 
ing problems due to losses from beam sputtering. The 
bombarding-cnergy-dependent correction factor F a (E p ) 
was minimal for our BE3 S detector data. Our preci- 
sion target profile measurements determined the effects 
of proton beam energy averaging in the target at each 
bombarding energy, and provided essential information 
on the stability of the target over time. 

We showed earlier (see Fig ll3|) that our present 
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FIG. 18: S-factor data from direct experiments, all normalized 
to a common value of 817(0) (the mean DB best-fit value of 
21.1 eV b - see Table EH). The error bars shown are relative, 
and do not include scale-factor uncertainties. Solid curve: DB 
plus a 1 + resonance with parameters determined from fitting 
our BE1 data. Dashed curve: DB only. Calculations and data 
are normalized in the energy range E crn < 1200 keV. 

7 Be(p,7) 8 B results agree well with DB plus a fitted 1 + 
resonance below 1200 keV or so. In Fig. ^| we compare 
the energy dependence of our results with the data from 
the other 4 modern direct experiments, all plotted with 
a common normalization corresponding to Si7(0) = 21.1 
eV b based on our DB (plus a 1 + resonance) fitted to data 
below 1200 keV (see Table EI3- From Fig. 1151 it can be 
seen that all direct experiments are in reasonable agree- 
ment on the energy dependence of the 7 Be(p,7) 8 B cross 
section in the energy range E cm < 1200 keV. This pro- 
vides an important basis for comparing the direct results 
with indirect Coulomb-dissociation experiments (see be- 
low). 

C. Comparison with indirect experiments 

1. Coulomb dissociation 

Coulomb dissociation (CD) experiments, in which a 
secondary radioactive beam of 8 B nuclei is dissociated 
into 7 Be + p in the field of a heavy nucleus such 
as Au or Pb. have now been performed by several 
groups 0, 0, 0, 0, |2(j. In these experiments, one 
attempts to avoid nuclear contributions by measuring 
at very small 8 B scattering angles. Corrections for the 
branch to 7 Be*(429 keV) are required. The relative 
weighting of El, Ml and E2 multipolarities in the virtual- 
photon spectrum responsible for the Coulomb breakup is 
very different than the weighting in the direct photon- 
emission spectrum. Hence, knowledge of the virtual- 
photon spectrum as well as the multipole decomposition 
of the direct 7 Be(p,7) 8 B cross section are both needed 



to translate the measured breakup cross section into an 
equivalent direct cross section. This is done with the as- 
sumption that the direct cross section, excluding the 1 + 
resonance, is pure El. In the virtual-photon spectrum, 
the E2 cross section, which is negligible in the direct pro- 
cess, is enhanced relative to El by several orders of mag- 
nitude and may not be negligible in the breakup cross sec- 
tion. E2 contributions estimated from measured breakup 
momentum distributions range from small but signifi- 
cant to negligible [2(j, and are not given reliably 
by theory. The effect of the 630 keV 1 + Ml resonance 
is taken into account, while weaker Ml strength located 
at higher energies (see e.g. the 3 + resonance discussion 
in Sec. lIVBl above) is not treated explicitly. Three-body 
Coulomb postdecay acceleration effects are a concern, es- 
pecially at low relative 7 Be + p energy and low bombard- 
ing energy. A recent estimate |4g suggests this effect may 
not be important in the work of refs. 0,H3, while for 
the 81 MeV/nucleon, small-angle data of ref. [l^, cor- 
rections for this effect have opposite signs for the points 
at w 200 and 400 keV relative 7 Be + p energy. 

The energy dependence of the inferred El S'17-factors 
for the 7 Be(p,7) 8 B reaction are shown in the top panel 
of Fig. ^3 determined from the 4 most recent CD ex- 
periments. To display the energy dependence clearly, we 
normalized each CD experiment to a common value of 
5i7(0) determined by fitting each CD data set below 400 
keV with the DB theory (see Fig. 4a of 20] for a sim- 
ilar plot with the original CD data, not renormalized) . 
The DB and Typel theory curves are also shown, again 
normalized to the CD data below 400 keV. As can be 
seen, the energy dependence of the inferred 7 Be(p,7) 8 B 
cross sections is similar in the works of Iwasa, Davids 
and Schumann, and is similar to the Typel theory, while 
the results of Kikuchi show a somewhat flatter energy 
dependence. 

It should be noted that theory cannot predict reliably 
the energy dependence of the 7 Be(p,7) 8 B cross section. 
Hence the Typel calculation is useful here only as a guide 
to the eye that reproduces the energy dependence of the 
CD data. 

We showed above that our direct measurements be- 
low 1200 keV have an energy dependence that is well- 
described by the DB theory (plus a 1 + resonance) except 
at the highest energies (1100 - 1200 keV) where DB falls 
a few percent low. The energy dependence inferred from 
the CD experiments is significantly steeper than the DB 
theory and does not agree with the direct results. The 
bottom panel of Fig. ^3 shows the CD data, not renor- 
malized, along with our present direct results (with the 
1 + resonance subtracted) and the best-fit DB curve to 
the present data. 

Because of the different energy dependences observed 
in CD and direct experiments, it is difficult to know how 
to proceed with a more quantitative comparison. If we 
ignore this problem, and focus on CD data below 400 
keV, as we did with the direct experiments, then our 
DB fits yield the 6*17(0) values shown in Fig. 1201 These 
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FIG. 19: El 7 Be(p, 7 ) 8 B S-factors inferred from Coulomb dis- 
sociation (CD) experiments. Bottom panel: absolute CD S- 
factors, together with our direct results (with the 1 + reso- 
nance subtracted) and the best-fit DB curve to our direct 
low-energy data. Top panel: CD data plotted with a com- 
mon normalization based on the mean value of 19.3 eV b for 
Si7(0) determined by fitting each data set to the DB theory 
below 400 keV. Solid curve: DB calculation; dashed curve: 
Typel calculation. The experimental error bars shown in all 
cases are relative, and do not include scale-factor uncertain- 
ties. 



values are mutually consistent, with a mean of 19. 3± 0.7 
eV b. A fit to this mean value together with the mean 
value deduced from direct experiments of 21.4 ± 0.5 eV 
b has probability P ~ 0.02 that these results arise from 
the same parent distribution. On the other hand, if we 
fit the CD data between 750 - 1400 (or 1000 - 1200) keV 
with the DB theory, the mean CD value is 6*17(0) « 22 
eV b, in very good agreement with the direct result (see 
Fig. ^3 bottom panel). However, there seems to be no 
independent motivation for fitting only high-energy CD 
data. 
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FIG. 20: CD ^17(0) values from DB fits to inferred 
7 Be(p,7) 8 B cross sections below 400 keV. The horizontal solid 
and dashed lines indicate the mean value 617(0) = 19.3 ± 0.7 
eV b. 



2. Heavy-ton transfer and breakup 

A Texas A&M group has used measurements of pe- 
ripheral heavy-ion transfer and breakup cross-sections to 
deduce the asymptotic normalization coefficient for the 
7 Be + p component of the 8 B ground state wavefunction. 
This coefficient, together with a capture- model calcula- 
tion (and an assumed P3/2/P1/2 ratio in 7 Be + p) can be 
used to infer Si 7 (0). The value Si 7 (0) = 17.3 ± 1.8 eV b 
has recently been inferred from the weighted average of 
10 B( 7 Be, 8 B)fBSe and 14 N( 7 Be, 8 B) 13 C results at E( 7 Be) 
= 85 MeV |2l|, and a variety of peripheral heavy ion 
breakup results at 28 to 285 MeV/u have been used to 
infer Si 7 (0) = 17.4 ± 1.5 eV b H3- However, a different 
analysis |47j of the same breakup reaction measured with 
a C target leads to a substantially larger Si 7 (0) value of 
21.2 ±1.3 eV b in good agreement with the direct mean 
value of 21.4 eV b. A determination of the asymptotic 
normalization coefficients for the P1/2 and P3/2 compo- 
nents of 8 Li — > 7 Li + n together with the assumption of 
mirror symmetry leads to Si 7 (0) = 17.6 ± 1.7 eV b for 
7 Bc + p |4^. These S- factor determinations thus tend 
to be even smaller than those deduced from CD experi- 
ments. 

S(E) values inferred from 16 0( 3 He,d) 17 F cross section 
measurements have been compared to direct 16 0(p, 
7) 17 F cross section measurements HilEI for capture to 
the ground state and to the first excited state. For the 
first-excited state transition, the (p, 7) results of ref. 
and the transfer-reaction results agree to approximately 
(6 ± 11)%, where the uncertainty is determined by the ± 
10% systematic uncertainty in the transfer reaction and 
± 5% uncertainty in the absolute (p, 7) cross section. For 
the ground-state transition, the central values from (p, 
7) |51jand from ( 3 He,d) agree within 10% or so but it is 
difficult to quantify the significance of the comparison |49| 
since the absolute (p,7)cross section uncertainty was not 
specified. 
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VI. RECOMMENDED VALUES FOR S 17 (0) AND 
Si 7 (20) 

In comparing results from experiments that employ 
very different techniques, it is necessary to group the re- 
sults according to technique, as we have done above, to 
see if the results are technique-dependent. The evidence 
presented here points strongly toward such a technique 
dependence. 

One could obtain a "best" or recommended value for 
Si7(0) by combining the mean results from different tech- 
niques, expanding the error on each mean by a common 
factor such that the combined fit has probability P = 
50%; i.e. x 2 /^ = 0.46 for v = 1. With this procedure, the 
combination of direct plus low-energy CD results yields 
Si 7 (0) = 20.7 ± 1.5 (expt) ± 0.6 (theor) eV b. Including 
the heavy ion reaction results as a separate, third group 
would increase the overall uncertainty and further lower 
the central value. This procedure has the disadvantage 
that it treats the results from different techniques on an 
equal footing. We have shown that the 7 Be(p,7) 8 B en- 
ergy dependence inferred from CD experiments disagrees 
with direct results. In our opinion, this difference in en- 
ergy dependence must be understood before results from 
direct and indirect experiments can be combined. In ad- 
dition, considerable theoretical modeling is necessary to 
infer the 7 Be(p,7) 8 B cross section from CD experiments, 
or Si7(0) from heavy-ion transfer and breakup experi- 
ments, and it is difficult to understand all the uncertain- 
ties associated with this modeling. 

It is important to note that neither indirect technique 
has been tested by comparison to a known direct result 
with sufficient precision to demonstrate that systematic 
uncertainties are understood at the level of ± 3 - 5%. 
We conclude that, at present, the indirect experiments 
are not sufficiently understood to be included in the de- 
termination of a recommended value. 

A new direct 7 Be(p,7) 8 B cross section measurement 
using a 7 Be beam would be useful as an independent de- 
termination of the absolute cross section with systematic 
uncertainties different from those of 7 Be target experi- 
ments. However, to make a significant contribution, such 
a measurement would have to have a total experimental 
uncertainty of 5% or better. 

We base our recommendation for Si7(0) on the mean 
of direct experiments as given in Eg. 1141 above, which we 
repeat here: 

Si 7 (0) = 21.4 ± 0.5(expt) ± 0.6(theor) eV b, (15) 

where the quoted errors are la. Although our discus- 
sion focusses on Si 7 (0), it seems clear that theoretical 
extrapolation uncertainties should ultimately be reduced 
by using instead the S-factor at 20 keV 01 > close to the 
Gamow peak. Our recommendation for this quantity and 
its 1(7 uncertainty is 

Si 7 (20) = 20.6 ± 0.5(expt) ± 0.6(theor) eV b. (16) 



VII. SUMMARY 

We made new Be(p,7) 8 B cross section measurements 
that extend our earlier results to lower energy and have 
improved systematic errors. Based on our new data with 
E cm — 116 to 362 keV and the cluster-model theory of 
Descouvemont and Baye, we determine Si7(0) = 22.1 ± 
0.6 (expt) ± 0.6 (theor) eV b, where the theoretical (ex- 
trapolation) error is given by the rms deviation of Si7(0) 
values from 12 different theories fitted to the same data. 
Our new result is in excellent agreement with o ur p revi- 
ously published value of 22.3 ± 0.7(expt) eV b [14j . and 
supercedes it. 

We have fitted all published direct 7 Be(p,7) 8 B mea- 
surements at E cm < 425 keV and E cra < 1200 keV with 
the DB theory. The modern experiments (Filippone, 
Hammache, Strieder, Baby and this work) give very con- 
sistent results. For E cm < 425 keV , the combined fit 
yields Si 7 (0) = 21.4 ± 0.5 (expt) ± 0.6 (theor) eV b, 
with y^/v = 1.2 (P = 30%). At present, the uncertain- 
ties in Si7(0) from the experimental cross sections and 
from the theoretical extrapolation are nearly the same. 
This points to the importance of further theoretical work 
to reduce the extrapolation uncertainty. 

We have also examined the 4 recent Coulomb disso- 
ciation experiments. These experiments infer a steeper 
energy dependence for the 7 Be(p,7) 8 B cross section be- 
low 1200 - 1500 keV than is observed in the direct ex- 
periments. Since the energy dependence of the present 
direct measurements is unambiguous (at the level of a 
few percent), this indicates a systematic error in the in- 
terpretation of the CD experiments. Fitting the CD data 
below 400 keV with the DB theory leads to a mean value 
of S i7 (0) = 19.3 ± 0.7 (expt) ± 0.6 (theor) eV b, which 
is only compatible with the mean of direct measurements 
at the level of 2%. Fitting only the high-energy region 
of these data yields Si7(0) ~ 22 eV b, in excellent agree- 
ment with the direct result; however, it is not clear that 
such a restricted fit is well-motivated. Peripheral heavy- 
ion transfer and breakup experiments lead to Si7(0) « 
17.4 eV b. 

Our recommended values for Si 7 (0) and Si 7 (20), based 
on the totality of existing data, are given in Eas. 1151 and 

m 



VIII. IMPACT OF NEW Si 7 (0) ON THE SOLAR 
MODEL AND ON NEUTRINO PHYSICS 

Our recommended Si7(0) has a combined experimental 
plus theoretical uncertainty of ± 0.8 eV b, or ± 4%. This 
represents a considerable improvement over the 1998 rec- 
ommendation by Adelberger et al. @] of Si 7 (0) = 19 ±2 
eV b (lc). It also represents a considerable improvement 
compared to Si 7 (0) = 19 ±%jl eV b assumed in BP00 [53. 
If incorporated into the BP00 solar model calculation [3] , 
the uncertainty in Si7(0) would no longer make an impor- 
tant contribution to the overall uncertainty in the calcu- 
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lated solar neutrino production rate from 8 B decay, and 
the overall uncertainty would be reduced from ± 17% to 
± 14%. 

Recent combined analyses of solar-neutrino plus ini- 
tial KamLAND results have limited the allowed oscil- 
lation parameters to the LMA (large mixing angle) re- 
gion [UEj]. These analyses, which are independent of 
the solar model, determine /b, total = 1-00 ± 0.06 [Hif . 
or 1.05 -1.08 [53j (no uncertainty quoted), depending on 
method of analysis, where /b, total is the total 8 B neu- 
trino flux (at the surface of the earth) in units of the 
BP00 flux. Since the SSM 8 B neutrino flux depends lin- 
early on Si7(0), incorporating our recommended Su(0) 
value and uncertainty into the SSM leads to a new SSM 
8 B neutrino flux /b, total = 1.13 ± 0.16, in units of the 
BP00 flux. Thus the correctness of the SSM has been 
dramatically confirmed. 

Within 3(7, these solar-model-independent analyses 
also restrict the allowed LMA region to a primary mini- 
mum at Am 2 w 7 x 10 -5 eV 2 and a secondary minimum 
at Am 2 w 1.5 x 10~ 4 eV 2 . In the analysis of ref. 
these minima lie at /b, total = 1-00 and 0.88, respectively 
(see Table 3 and Fig. 4 of 54]). Thus the disfavored 
(secondary) minimum with /b, total = 0.88 is somewhat 
disfavored additionally by the SSM rate calculated with 
our recommended Si7(0). 

Because the combined atmospheric neutrino, solar 



neutrino-plus-KamLAND, and LSND results require 3 
distinct regions of allowed Am 2 , a fourth (sterile) light 
neutrino is required if all 3 data sets are correct. How- 
ever, even if the LSND results turn out to be spurious, 
sterile neutrinos may exist and play important roles in 
nature (see e.g. Existing limits on sterile neu- 

trinos are derived from solar-neutrino plus Kamland- 
antineutrino data |54j assuming CPT invariance. But the 
CPT properties of sterile neutrinos, unlike those of other 
particles, have not been tested. A rigorous test for sterile 
neutrinos requires knowing the production and detection 
rates for at least two neutrino sources with substantially 
different energy spectra. These could be the pp neutri- 
nos (whose production rate is determined by the solar 
luminosity) and the solar 8 B neutrinos whose production 
rate is affected by the results of this work. The new SSM 
8 B flux based on our recommended Si7(0) would be use- 
ful in an analysis that tests for sterile neutrinos without 
assuming CPT invariance. 
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